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We study the impact of flavour in thermal leptogenesis, including the quantum oscillations 
I of the asymmetries in lepton flavour space. In the Boltzmann equations we find different 

numerical factors and additional terms which can affect the results significantly. The 
upper bound on the CP asymmetry in a specific flavour is weaker than the bound on 
the sum. This suggests that - when flavour dynamics is included - there is no model- 
independent limit on the light neutrino mass scale, and that the lower bound on the reheat 
temperature is relaxed by a factor ~ (3 — 10). 



o 



Oh 

CD 



X 



PACS number(s): 12.60. Jv, 14.60.Pq, 11.30.Fs 

Leptogenesis, Neutrino Physics, Baryon Asymmetry of the Universe 

Note added: The equations (|lip and (|19|) of the published paper, which include the decaying 
oscillations in flavour space, do not reflect the flavour structure derived in the Appendix, or reproduce 
the behaviour discussed in the text. In this version, we modify equations (fTTI) and (fl9j) . and add a 
discussion, in Appendix A, on how to obtain these improved equations. We thank Georg Raffelt for 
drawing our attention to this, (and also to the paper by N.F.Bell et al, Phys. Lett. B500 (2001) 16, 
which indicates that fast gauge interactions should not affect the coherence of flavour oscillations). 
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the change of form of equations (jlip and (|19p . However, (|2QI2ip are are wrong by the factor z. 
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1 Introduction 



Leptogenesis pQ is an attractive mechanism to explain the baryon asymmetry of the Universe (BAU) 
[2]. It consists of the dynamical production of a lepton asymmetry in the context of a given model, 
which then can be converted to a baryon asymmetry due to (B + L)-violating sphaleron interactions 
|3j which exist in the Standard Model (SM). 

A simple model in which this mechanism can be implemented is "Seesaw" (type I) [4], consisting 
of the Standard Model (SM) plus 3 right-handed (RH) Majorana neutrinos. In this simple extension 
of the SM, the usual scenario that is explored consists of a hierarchical spectrum for the RH neutrinos, 
such that the lightest of the RH neutrinos is produced by thermal scattering after inflation, and decays 
out-of-equilibrium in a lepton number and CP-violating way, thus satisfying Sakharov's constraints. 
In recent years, a lot of work [51 [6], has been devoted to a thorough analysis of this model, the 
constraints that can be imposed, the relationships that arise with low energy neutrino physics, etc. 
For instance, a lower bound on the reheating temperature Trh and an upper bound on the scale of 
light neutrino masses can be obtained [71 [HI [9| [TO] . 

However, in order to obtain an accurate value of the produced lepton asymmetry it is necessary 
to solve (in a first approach at least) the system of coupled Boltzmann equations (BE) which relate 
the abundance of the lightest RH neutrino with the asymmetry produced in lepton number. The 
effects of the charged lepton Yukawa coupling matrix [h] are usually not considered. In this study, 
we take into account these effects in deriving the evolution equations for the asymmetries in each 
flavour, and study the implications for leptogenesis. Flavour effects have been also discussed in 
More recently, the generation of the family lepton asymmetries has been discussed in [12] and |13j . 
where the case of leptogenesis in the minimal seesaw model with two heavy Majorana neutrinos was 
considered; the role of flavours in the leptogenesis scenario from N2 decay in discussed in [14] and 
for "resonant leptogenesis" (degenerate Nj) in [15]. In this work, we find three effects: a qualitative 
difference in the equations, which allows quantum oscillations of the asymmetry in flavour space. 
We naively anticipate this could change the final asymmetry by a factor ~ 2 |16| . Secondly, we 
find different numerical factors compared to the case when no flavour issues are included in the 
Boltzmann equations, and additional terms when the flavour-dependent decay rates and asymmetries 
are unequal. Finally, we find that the CP asymmetries in individual flavours can be larger than the 
sum, which affects the leptogenesis bounds on Trh and makes the bound on the light neutrino mass 
scale disappear. 

In section [2] we review our notation and summarise the usual way the calculation is performed. In 
section[3j we present equations for the lepton asymmetry, represented as a matrix in flavour space: the 
diagonal elements are the flavour asymmetries, and the off-diagonals encode quantum correlations. 
These could be relevant, because, as will be explained, the muon Yukawa coupling comes into thermal 
equilibrium around T ~ 10 9 GeV, allowing the possibility that the asymmetry "oscillates" in flavour 
space. To motivate the equations used in section [3] we study a quantum mechanical toy model of 
simple harmonic oscillators. This study is explained in appendix A of this note. We concentrate 
on the BE for the three flavour asymmetries in section HJ and compare the resulting equations for 
the lepton asymmetry with the usual approximation. We present results in section [5l which aims to 
be comparatively self-contained. Appendix B contains a brief discussion on the AL = 2 terms, and 
Appendix C relates the asymmetries in 5/3 — Lj to those in the lepton doublets. 

2 Notation, Review and the "Standard" Calculation 

We consider the Lagrangian for the lepton sector 

C = h ik He Ri £ k + \ Jj H c Nj£j + ^N K N K + h.c. (1) 
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in the mass eigenstate bases of the singlet ("right-handed") particles N and e#. The mass eigenstate 
basis for the doublet leptons will be temperature dependent, so we do not fix it here. We assume 
from now on a hierarchical spectrum for the RH neutrinos, Mi <C M2 < M3. 

The light neutrino mass matrix, in the charged lepton mass eigenstate basis, is 

m v = U*D mv U ] = X T M~ 1 Xv 2 , (2) 

where U is the PMNS matrix, Da is a diagonal matrix of real eigenvalues of A and v the Higgs 
vacuum expectation value. 

The physical process of interest is the production of a lepton asymmetry in the decay of the 
lightest singlet neutrino N\, in the thermalised early Universe. We make the usual assumption, 
that interactions that are much faster than the expansion rate H — such as gauge interactions, 
{B + L)-violation, and (possibly) the r- Yukawa coupling — are in equilibrium, which imposes certain 
conditions [17] on the distribution functions of particles. For instance [18} I19j. when 

r T ~ 5 x 10~ 3 /£T > H (3) 

the chemical potentials of the singlet and doublet rs satisfy the relation + \i TL = fi TR . We neglect 
interactions whose rates are much smaller than H, e.g. those induced by the electron Yukawa coupling. 
Interactions that are of order the expansion rate, and those that are lepton number violating, are 
included in the evolution equations for the number density (eventually, Boltzmann equations). These 
are the iVi-interactions, the AL = 1 and AL = 2 scatterings, and possibly the interactions with 
strength set by the fx- Yukawa coupling, h u . 

It is straightforward to write down Boltzmann equations for the conserved asymmetries of the 
Standard Model (SM). The SM interactions which are in equilibrium define the conserved quantum 
numbers. For instance, below T ~ 10 5 GeV, when all interactions are in equilibrium, the conserved 
quantum numbers are the {B/3 — L,} for % = e,fj,,r. The Boltzmann equations are obtained by 
identifying processes that change the asymmetries, and summing their rates. The only delicate point 
arises in the temperature range where the interactions induced by tin come into equilibrium, because 
this changes the relation between and • 

The Boltzmann equations for number densities neglect the possibility of quantum effects, such as 
oscillations, due to interference among different processes. Such effects could appear in the equations 
of motion for the (flavour-dependent) number operator 

flip) = 4(p)a{(p) - a£(p)al(jt), (4) 

which counts the asymmetry in lepton doublets^. Notice the inverted flavour order between the 
particle and antiparticle number operators [20]. The operator is a matrix in flavour space, analogous 
to the density matrix of quantum mechanics, and is sometimes called a density matrix. When the 
flavour indices ij are the charged lepton mass eigenstates, then the diagonal elements are the flavour 
asymmetries stored in the lepton doublets, and should satisfy Boltzmann(-like) Equations. The 
trace, which is flavour-basis-independent, is the total lepton asymmetry. The off-diagonals encode 
(quantum) correlations between the different flavour asymmetries, and should decay away when the 
charged Yukawa couplings are in equilibrium [llj. Variants on such an operator have been studied 
in the context of neutrino oscillations in the early Universe [20\ l2Tj \22\ I23j . and in particular the 
generation of a lepton asymmetry by active sterile oscillations [24\ I25j . This operator makes brief 
appearances in [11] where the issue of flavour in the context of thermal leptogenesis was first discussed. 

1 / wears a hat as an operator, to distinguish it from its expectation value which we will usually discuss in the 
Appendix. The operator a+ (a_) creates particles (anti-particles) of flavour i. 
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In appendix we will motivate equations of motion for these flavour-dependent lepton number 
operators using quantum mechanical oscillators. For simplicity, we separate the (possibly quantum) 
flavour effects from the (assumed classical) particle dynamics. That is, we extract the flavour structure 
from the toy model, and input the particle dynamics and the Universe expansion, by analogy with 
the Boltzmann equations. This is an elegant way to obtain flavour-dependent Boltzmann equations, 
because the tensor structure of the operators in flavour space is clear. In a later publication, we 
will derive equations of motion for /V in field theory. Indeed, its dynamics could be obtained from 
the transport equations deduced from the finite temperature Schwinger-Dyson equations involving 
flavoured fermionic propagators. In this set-up [26], one makes use of the closed time-path formalism 
to describe nonequilibrium phenomena in field theory leading to a complete nonequilibrium quantum 
kinetic theory. The operators f£ are associated to propagators which may mix the flavour and which 
are "sourced" in the Schwinger-Dyson equations by the flavour-dependent CP asymmetry €ij from 
the decay of the right-handed neutrino. 

The usual starting point for leptogenesis is a system of coupled equations. If the asymmetry 
is due to the decay of the lightest RH neutrino Ni, and the processes included are decays, inverse 
decays and lepton number violating scattering (AL = 1 from Higgs exchange, and AL = 2 due to 
RH neutrinos), then these equations are: 



dY L 

dz 



dY Nl 
dz 



sH(Mx 



sH(M 1 



{•JD + 7AL= 



Y Nl 

Y N eq 



1 



Yn, 



\ Y L . 

1 I ejD — yeq KID + 7AL=1 + 7AL= 



(5) 



(6) 



IjVi is the abundance of the RH neutrino N\ normalized to the entropy density s, and Yl is the 
sum over flavour of the difference between normalized abundances of leptons and antileptons. In 
eqs. © and ©, Y? q is the equilibrium number density of a particle i, z = -yr- and the quantities 
7r>! 7AL=i and 7al=2 are thermally averaged rates of decays, inverse decays, and scattering processes 
with AL = 1,2 and include all contributions summed over flavour (s, t channel interference etc). The 
explicit expressions are given in the literature (see for example [5]) and we will discuss them when 
needed further on. The total CP-asymmetry e produced in the decay of the lightest RH neutrino is 
given by [27] 

' ^ * ' ' ' (7) 



8vr[AAt] 

where / is the loop factor [27J, given by 

m = 



ii 



[\tf]ljf(Ml/M?), 



,-> \ i 1 + x 
(1 + x) In h 



(8) 



We recall that the rates d and 7al=i are both proportional to the tree-level decay rate given by 
and the out of equilibrium condition is satisfied when: 



T=Mi 



^|<l. 



(10) 



We exploit this expression in the next section. 
3 For a later use, we give the expression of the reaction density jd = kI(z\ ^' wnere K n (z) are the K type Bessel 
functions and is the equilibrium density. 
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In other words, comparing fh\ torn* ~ 2.3 xlO" 3 eV gives a measure of whether the decay is occurring 
out of equilibrium or not. The rescaled decay rate rh\ is bounded below: m\ < rh\, and naturally of 
order y Am^ tm when the light neutrino spectrum is hierarchical [28j . 

The term 7al=£I has a more complicated dependence on the masses of the neutrinos, and can be 
treated in two separate regimes. The contribution in the temperature range from z = 1 to z ~ 10 is 
proportional to rh\ while for the range of temperature (z> 1), the dominant contribution at leading 
order is proportional to to 2 = Tr(mtwv) = Si m i- Here denotes the mass eigenvalues of the 
left-handed neutrinos. 



3 Flavour Oscillations 

In this section we consider thermal leptogenesis for two flavours, which to be explicit, we label to 
be the e- and the \x leptons. This is interesting in two cases. First, if leptogenesis takes place when 
<C H (see eq. p])), it will be described by Boltzmann equations for the r and number densities, 
where i\ = l\ — (£\ ■ f)f and = Xu is the decay direction, in lepton flavour space, of N%. So the 
equations of this section could be applied in this high temperature case, with the relabelling \x — > r, 
e — ► £'i. Secondly, if the muon Yukawa coupling comes into equilibrium as the asymmetry is being 
created, there could be oscillations of the asymmetry, parametrised in off-diagonal elements of the 
asymmetry density. Focusing on this quantum effect, we neglect the third lepton generation, because 
we suppose we are in the interesting range of temperatures ^ 10 12 GeV where the interactions induced 
by the r- Yukawa coupling are in chemical equilibrium. Therefore no quantum correlations between 
the lepton asymmetries involving the third family are expected to survive. Adding the third family 
will be discussed in the next section. 

The equations describing the two-flavour system have been obtained in appendix [A] to which 
we refer the reader for more details. To clearly see the effect on flavour oscillations due to the 
muon Yukawa coupling coming into equilibrium, we map eq. (|46p onto two-flavour equations for the 
asymmetry in the early Universe, neglecting the transformation to (B/3 — Li). As discussed after 
eq. (fl6|) . these equations describe decays and inverse decays of N\, and the resonant part of AL = 2 
scattering is also included. We work in the charged lepton mass eigenstate basis, and distribute 
indices by analogy with the model of the appendix (f%i «-> Y^ , and A J A* J <-> 7 U ). 

We assume, to obtain these equations, that the flavour-blind gauge interactions can be neglected. 
However, the oscillation frequency in flavour space depends on the energy of the lepton, and within 
the oscillation timescale, a lepton will participate in many energy-changing gauge interactional). In 
this note, we use the thermally averaged energy (E) to estimate the oscillation frequency |24j . That 
is, we approximate the integral (i J Edr) along the path from one lepton-number violating interaction 
to the next, to be i{E) f dr. We will include the gauge interactions more correctly in a later analysis 

The system describing the two flavours is given by 



dz 



yee 
L L 



sH{M{) 



ID 



yeq 



2Y[ q 



Id 
^d 



I'd" 

HZ 



ye, 
ype y^n 



IK 



Vee 

L 
V fj,e 



ye, 
1 L 



mi 



Vee 



ye. 



4 We emphasize that — g7 D . corresponding to RIS (real intermediate states) in the AL = 2 interactions have to be 
carefully subtracted, to avoid double counting in the Boltzmann equations [5]- 
5 We thank A. Strumia for discussions about this point. 
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where the parenthesis [-, •] and {•, •} stand for commutators and anti-commutators, respectively. The 
matrix [773] is defined as 



lk\ 



(12) 



where 773 = YlilD * s the total (thermally averaged) decay rate. For i = j, is the decay rate of 
Ni to the i-th flavour. The CP asymmetry in the i-th flavour is e M , e = Y^i £M > an d notice that e %3 is 
normalised by the total decay rate (so that [e] transforms as a tensor under i basis rotations). The 
matrix [e] is now defined as 



v 3 



(16vr) [AAt] n y 



{(A lj )[AAt] 1J A} J -(A^.)[A*A T ] 1J A Jl } 



J 1 M 2 



Mi 



(16th; 2 ) [AAt] n 

where f(x) is the loop function eq. 
The matrix A is given by 



Im{AiA CT [m*] CT j - A* [m u ] ai X*} 



My 
1 

(for Mi < M 2 ,M 3 ) 



(13) 



A 



(14) 



A ee 
A«" 

where A ee and A w indicate the complex thermal masses due to the electron and muon Yukawa 
couplings: 



These will cause the decaying flavour oscillations of the asymmetries. Since the interactions mediated 
by the electron Yukawa coupling are out-of-equilibrium, we may safely set A ee = 0. On the other 
hand [TS] 



h 2 



fifi — T, ^fifi — 5 x 10 hp T, 



(16) 



leading to 



Af ^3xHT 4 /^, 
d » Mi 



(17) 



where Mp is the Planck mass. Interactions are mediated by the muon Yukawa coupling only for 
z ^ Zd ~ 1/A^\ Therefore, as thermal leptogenesis takes place at temperatures on the order of Mi, 
we conclude that - roughly - interactions involving the muon Yukawa coupling are out-of-equilibrium 
in the primeval plasma if Mi <; 10 9 GeV. 

Following Stodolsky [21], we may now write eq. (jlip in a more compact form by expanding the 
various matrices on the basis provided by the <T-matrices, = (I, a). A generic 2x2 matrix A can 
be written as 

A = A^ = \TT{Aa^}. (18) 

In this basis, the Boltzmann equation (jlip becomes a system of equations of the form 



dz 

dY L 

dz 



Y, 



sH(Mi 



sH{M 1 



ID 



ID 



1 



l D Yl 



1 

1 L 



ID ■ Y L 



v e i I 

, Ni ) 

- 1 J e- jD^g ~ ^IdYl) + 4 x Y L + A d x A d x Y L . (19) 
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Notice, in particular, that Y® = (l/2)Tr Yj, and therefore it represents half of the total lepton asym- 
metry. On the other hand, Y£ = (1/2) (Y ee — Y^) indicates half of the difference of the asymmetries 
in the electron and muon flavour densities. The components Y£ and YV parametrize the off-diagonal 
entries of the asymmetries and therefore account for the quantum correlations between the lepton 
asymmetries. As the matrix A in eq. (|14p is diagonal, only the A and A z components are different 
from zero. 

A simple inspection of eqs. (|19p tells us that the components Y£ and Yf precess around the z- 
direction with an angular velocity set by the thermal mass to^. At the same time, such a precession 
is damped by the interactions mediated by the muon- Yukawa coupling at a rate ~ L„. Notice 
that ujfj.fi ~ r w , so unlike the case of neutrino oscillations in matter (MSW), the decoherence and 
oscillation timescales are the same. The term A x Yl contains all the information about the action 
of the decoherent plasma onto the coherence of the flavour oscillations: if the damping rate r w is 
much larger than the expansion rate of the Universe, the quantum correlations among the flavours 
asymmetries are quickly damped away. Therefore, we expect that, if leptogenesis takes place at values 
of z 3> 2d) quantum correlations play no role in the dynamics of leptogenesis [11] . 

If thermal leptogenesis occurs well before the muon- Yukawa coupling enters into equilbrium (that 
is at z <C Zd , T ~ Mi 3> 10 9 GeV), then \A\z <C 1 and the last term of the second equation of (fT9l) can 
be safely dropped. In the absence of the charged lepton Yukawa coupling, we may rotate the system 
in such a way to put all the asymmetry generated by the decay of the right-handed neutrino into a 
single flavour, e.g. Y ee . Under these circumstances, Y® = Y£ = (Y ee /2), Y£ = Y^ = and similarly 
for the CP-asymmetries, e° = e 2 = (e ee /2), e x = e y = 0. By simply summing the corresponding 
euqations for Y® and Y£, one immediately finds the Boltzmann equation for the single flavour case. 

To get some insight about the dynamics of the system, let us now discuss the solution of eqs. (|19p 
in the case in which thermal leptogenesis takes place nearly in equilibrium which, given the low-energy 
neutrino parameters, is likely to be the one realized in Nature. We parametrise the inverse decays with 
rates 7^ = n£ (K x {z) / K 2 (z)) T ij through the parameters (defined in eq. (pH) K ij = (T ij /H) T=Mi 
and work in the limit in which every ^ 1. By defining the combinations Y^ = (Y£ ± iYj\ /y/2 
and by exploiting various saddle-point approximations, we find 



Yl * 






9* 




e z 


Yl * 






g* 




e± 


y£ * 






9* 



Vf ^ - Jf6 ^1 + K+Y+ + K-Y L ) («), 



2K°z { K° L[ f) > 
;T 3/2 n - z& TA,T 



z d e Q cos 

,5/2 



2 



2 2\ 



K^Y^z^e^cos 



2 

r aw 



e 2 d 



2 2\ 

z f - Zd) 



M7{4-*1) } (20) 



where Zf ~ \n.K° + {h/2)\n.\n.K° . These solutions have been obtained in the limit in which z& & Zf, 
that is when the muon- Yukawa coupling induces interactions that are in thermal equilibrium at 
temperatures larger than the effective temperature at which inverse decays go out-of-equilibrium. In 
the limit z^ <C zi , eqs. (f20l) clearly show the decoupling of the quantum correlations from the dynamics 
in the limit in which muon Yukawa coupling induced interactions have come into equilibrium long 
before thermal leptogenesis, A^Zf S> 1. In this limit, we find 
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As a last comment, let us notice that in the limit Zf — z^, the quantum correlations are not efficiently 
damped out, and appear as oscillations in the washout terms of eq. (fl"9j) . We expect that the effect 
on the final total lepton asymmetry, 

Y? oc — , 22 

(K°) 2 -\K\ 2 

will be analogous to multiplying the washout rate by a neutrino survival probability, and integrating 
over time (an 0(1) effect). This case deserves further investigation and a detailed numerical study 



4 Flavour in the Boltzmann Equations 

In this section, we consider Boltzmann equations (BE) for the diagonal elements of [Yl] with 3 
flavours. As discussed in the previous section and in [TT], BE are appropriate when the interactions 
of the charged lepton Yukawa couplings are much faster, or much slower, than the expansion H— 
provided one works in the physical basis. To be explicit, we work in the charged lepton mass eigenstate 
basis, so these equations apply when the r and /j, Yukawa couplings are in equilibrium (see eq. ([3])), 
and have caused the off-diagonal elements of [Yl] to be irrelevant. This corresponds to T ^ 10 9 GeV 
for the SM (and to T ^ 10 9 x tan 2 (3, where tan /3 is the ratio of the vacuum expectation values of 
the two Higgs bosons, in the supersymmetric version of the SM). 

Let us start out by considering the Boltzmann equation as a matrix equation in flavour space, with 
only decays and inverse decays, as we did in the previous section, and neglecting H the transformation 
to B/3-Li 



sH{M{) d_ 

z dz 



\ree 



ID 



2Y e L q 



Id 



-ft 



TT 



Yf 



VTT 



'(23) 



where e tl and 7^ are defined as in eqs. (|12l I13|) . Since the off-diagonal elements of the matrices have 
been neglected, this equation is not invariant under flavour rotations, and only applies in the physical 
basis where the off-diagonals are suppressed. 

We can illustrate the problem that arises when one sums over flavour by taking the trace of eq. 
pi . The BE for the sum is 



^ dz 



sH 



Y Nl 



1 ejD 



v 



,1D 



V e, 3 



(24) 



where e n = e is the total CP asymmetry. Like jd it is the trace of a matrix in flavour space, see 
eq. ([JJ). 

In the "single flavour or dominant state approximation" (see eq. ([6])), the trace of the product of 
matrices in the last term on the RHS of eq. (I24p is approximated as the product of the traces: this 
term is replaced by ■ Y[ J ^ Yd- This is correct in the absence of charged lepton Yukawa couplings; 



6 Including this transformation would mildly reduce the washout term, because, e.g., an £^ produced is a unit of 
B/3 — L M , but not all the B/3 — L M is in ^ M s. In appendix C we give the explicit tranformation from Li to B/3 — Li 
for different cases. 
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if I is the flavour combination into which N\ decays, then there is a basis where [jd] = diag(7£), 0, 0) 
and [Yl]ii = Yl. The trace is invariant under basis transformations in flavour space, so for any 
unitary V 

^Tr{V [7r>] V f , V[Y L ]V j } = j D Y L (no charged lepton Yukawas). (25) 

In particular, this is the washout term evaluated in the charged lepton mass eigenstate basis, see eq. 
©• 

If the charged lepton Yukawa couplings are included, then new terms appear in the equations for 
the asymmetry [II], as discussed in the previous section. These new terms suppress the off-diagonals 
of [Yi] in the charged lepton mass eigenstate basis, so the Boltzmann equations (without the off- 
diagonals) are only appropriate in the charged lepton mass eigenstate basis, and in general it is no 
longer true that Tr{[-f D ][Y L }) = Tr[j D ] Tr[Y L ]. 

One way to see what is being neglected, if the single flavour approximation is used when the 
charged lepton Yukawa couplings are in equilibrium, is to expand the diagonal matrices [Yl], [e] and 
["Yd] ° n the identity matrix I and the diagonal SU(3) generators A3 = Diag{ 1,-1,0} and As = 
■^Diagjl, 1, —2}. For instance, 

[Y L ] = Y I + ysAs + y§As (26) 

where Y 3 = ±Tr(YA 3 ) is the flavour difference (Y ee - Y^) /2. Taking the trace of eq. ((23J), we obtain 
the BE for the total asymmetry Yj,: 

dY L : ((Y Nl \ (1 Y L Y 3 Y 8 \\ 

-^ = 7H\\Y^ i - l ) lDe - U 7D ^ + lD "Yf + 2lDS Y? ) ) ■ (27) 

In the absence of AL = 2 interactions, eq. (|23p consists of three decoupled equations, so the final 
asymmetry can be obtained easily by solving for each flavour, then adding the solutions. However, 
we can also obtain Boltzmann equations for the asymmetries Y3 and Yjj, by multiplying eq. (|23p by 
A3 (or Ag), and then taking the trace: 

dY 3 z ( (Y Nl \ fr/D Y 9 7D3 Y L 7 D3 Y 8 7D8 Y 3 \ \ 

" m[{y*- V 7De3 "U^ + ~^ + 7f^ + 7f^JJ' (8) 

^^8 _ JL((X^_-\\ (lD Y % , 7£>8 Y L , Idz Y3 7Q8 Y% \ \ (1?Q , 

dz sH \\y% J 1D€8 V 3 Y? + 3 Y? + ^3 Y? V3 Y? ) ) ' 1 j 

The factor of 1/3, multiplying 7_dYl in eq. (|27p is interesting. One can see from [29J that in many 
models, N\ decays about equally to all flavours, and produces about equal asymmetry in all flavours. 
One can therefore neglect the terms 7D3Y3 + 7_dsYs in eqs. (f2Tl) and ([29]) . Thus, leptogenesis is 
described by the usual equation, with the washout reduced by a factor of 1/3. This should not be 
surprising: 1/3 of the asymmetry is in each flavour, and for each flavour, the (inverse) decay rate 
is 1/3 of the total. Summing over flavour, 3 x 1/3 x 1/3 = 1/3. We will study leptogenesis in the 
context of these flavoured Boltzmann equations in a subsequent publication [16J. 



5 Implications of flavours 

In this section, we study the modifications due to flavour of the leptogenesis bounds on Trh and the 
light neutrino mass scale in. In the "single flavour approximation", an N\ -parameter- independent 
bound on the light neutrino mass scale, of 4= ^ 0.15 eV, can be obtained [91 [10] . We will argue below 

that this bound does not hold when flavour effects are included. The bound m u ^ 4y / 10 10 GeV/Te_/ / 
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eV [3D], from requiring AL = 2 processes to not wash out the asymmetry, would apply, but remains 
a function of the (unknown) temperature Tb-l at which leptogenesis takes place. 

The limit m u 5s 0.15 eV, can be understood to arise from the lower bound on the total decay 
rate mi < mi, and the upper bound on the total CP asymmetry [7J[8]: 

where m.3 is the largest neutrino mass, and mi the smallest. We assume the light neutrinos are 
degenerate, so |mi| ~ \m2\ ^ \m^\ = fh/y/3 — but the masses can have different Majorana phases. 
Leptogenesis takes place in the strong washout regime, due to the lower bound on the total decay 
rate. The final baryon asymmetry can be roughly approximated as t]b ~ I0~ 3 e/K oc Am 2 tm /fh 2 . 
As the light neutrino mass scale is increased, Mi and the temperature of leptogenesis must increase 
to compensate the Am^/m 2 suppression. However, this temperature is bounded from above, from 
the requirement of having the AL = 2 processes out of equilibrium when leptogenesis takes place: 

m 2 T 3 , 10T 2 



12t™ 4 ^ M P ' (31) 

so M\ ^ 10 10 (eV/m) 2 GeV. There is therefore an upper bound on the baryon asymmetry which 
scales as 1/m 4 , and with our rough estimates, one finds m/V3 & 0.1 eV. 

The individual flavour asymmetries, en do not satisfy the bound of eq. (f3"0"j) . This can be guessed 
from the "Jarlskog invariant" for leptogenesis which is the trace: 

Jl = 9{Tr[m i ,mt(A T A*)- 1 m i ,(A t A)- 1 m t J} 

][U*DlW T DfW*D m ^D- 2 WD mv U T ]A , (32) 




where the unitary matrix W transforms from the basis where m u is diagonal to the one where A^A 
is diagonal (e.g. in the D\ basis, m v = W* D mv W^). The total asymmetry e is proportional to this 
trace, which vanishes when the light neutrinos are exactly degenerate in magnitude. However, if the 
invariant is "cut open" on lepton flavour indices, one can show that the flavour asymmetries satisfy 

M < -^ r , (33) 

and only the sum vanishes. So if we can arrange — — e TT ~ v / 3A/im/(87r'u 2 ), and decay rates 
li) = s2 lDi Iff = c 2 lDi then in the approximation of strong washout and neglecting AL = 2 
processes, eq. (f2T|) gives the total lepton asymmetry to be 



1 3Mim* c 4 - s 4 
g* Zf 8ttv 2 c 2 s 2 



So we obtain that the baryon asymmetry is independent of the light neutrino mass scale, under the 
assumption that AL = 2 processes are negligeable. In passing [§, one can notice that Yl increases, 
when one of the partial decay widths decreases, provided that it remains in the strong washout regime: 
s 2 , c 2 3> rh*/rh. 



7 This invariant is only applicable when the right-handed neutrinos Ni are hierarchical. 

8 We will study the final asymmetry as a function of similar angles in a subsequent publication [16] 
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We first show that flavour asymmetries are constrained by eq. ()33j) . when light neutrinos are 
degenerate at mass m / v3, as opposed to the stronger bound on the total asymmetry of eq. (f3"0"j) . For 
hierarchical Nj, with exactly degenerate light neutrinos 

-Im{Ai(A • mji} = — — = < — — =— , (35) 



8vn; 2 |A| 2 " 8ttv 2 |A| 2 ~ 8vrv 2 



where there is no sum on i, (A)j = Aij, and m* = UDmU 1 = ^UU T 

We can verify that a large neutrino mass scale is consistent with leptogenesis by explicit con- 
struction of a model. It is convenient to use the Casas- Ibarra |32j parametrisation, in terms of mj, 



— 1/2 — 1/2 

Mj, U and the complex orthogonal matrix R = vD M \UD m . This ensures we obtain the correct 
low-energy parameters. We take a two-flavour model, and write R as a rotation through the complex 
angle 9 = (p + irj. The flavour asymmetries, for maximal atmospheric mixing and no phase in U, are 

VSMifh 1 Cjt{T> rr* c> tt i VSMifh sinhr/coshr/ 

^ = 8™ 2 E^P^F { lk ^ 1P " P} = ~ 8- 2 sinh 2 , + coslV° S2 ^ 
VZMifh 1 v^Mim sinhr/coshr/ 

£ ^ = "5 IP D TT * \0 ^{ R lkU Tk RlpU Tp \ = — 5 —2 -5— COS 2(^9 . (36) 

87TV z2a\z2j R ij U aj\ 87ru 2 smh 2 rj + cosh 2 r/ 

The flavour dependent decay rates 7 M/J and 7 rT are proportional to 

2 Mim 2 Mim , . 2 2 

A = 2v 2 y/3 = ^-^= (smh t? + cosh r ? + sm2v?) , 

2 M\fh 2 M\m 2 2 

r = 2zj 2 ~ = 2u2y^ r? + cosh 77 -sin 2^). (37) 

We see that we can choose r\ and (p as required to obtain eq. (|34f) with c 4 — s 4 ~ s 2 c 2 . We conclude 
that, if Mi £ 10 12 GeV and flavours effects are relevant, the absolute bound on the light neutrino 
masses is the one inferred in ref. [30J, m u 4y / 10 10 GeV/Te_ £ eV. Of course, if M\ P> 10 12 GeV, 
flavour effects are irrelevant and we are back to the single flavour case where the bound ^ 0.15 

n 

eV applieqj. 

The minimum required reheat temperature of the Universe, Trh, will depend on the flavour 
structure of the lepton asymmetry. For thermal leptogenesis to take place, an adequate number 
density of N\ must be produced by thermal scattering in the plasma, suggesting Trh ?Z Mi/10 [HI Ej - 
There is a lower bound on M±, from requiring a large enough CP asymmetry; the corresponding 
bound on Trh depends on the fine details of reheating, whose dynamics have been addressed in ref. 
[5j. The lower bound of [5], Trh ^ (2 — 3) x 10 9 GeV, was obtained with hierarchical light neutrinos 
and strong washout. We include flavour effects by rescaling this bound. 

In many popular seesaw models, the light neutrinos are hierarchical, leptogenesis takes place in 
the strong washout regime K > 1, and the asymmetries and decay rates are similar for all flavours: 
7r> ~ 7d an d gM ~ eJJ ■ I n this case, e is bounded by (|30|) . but the washout is reduced by the factor 
of 1/3 from equation (|27|) . so we estimate that 

T RH £ 6 x 10 s GeV . (38) 

Thermal leptogenesis could work with a lower reheat temperature if the washout rate in one 
flavour is suppressed. If K > 1, a thermal distribution of N\s will be created, but the asymmetry in 
flavour j is only washed out by 7^, so0 Yp ~ 10~ 3 e^ /K^ . For instance, in our "degenerate" two 



9 We thank E. Nardi, Y. Nir and E. Roulet for discussions about this point. 
10 This formula is valid in the strong washout regime, so requires — 7" K/^yo > 1. 
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generation model, we can tune the couplings lj so that Xi^/Xir ~ s is of order fh^/fh where m ~ 
eV, in a generous intepretation of the cosmological bound [33]. This suppresses e MAt with respect to 
eq. (|34p by s, so the final asymmetry Yl ~ Y^ scales as 1/s and we obtain 

Trh ^ 2s x 10 9 GeV ~ 10 8 GeV for tuned, degenerate light neutrinos. (39) 

This appears theoretically unmotivated, compared to other leptogenesis scenarios at low Trh [T51I54"] . 

Our results can therefore relax the tension between thermal leptogenesis and the gravitino bound 
|35j in super symmetric extensions of the SM which impose an upper bound on Trh of about 10 9 GeV. 

6 Summary 

Including flavour in the analysis of leptogenesis modifies the equations of motion (Boltzmann Equa- 
tions) for the asymmetry, and changes some of the relations between parameters of the calculation. 

The interactions related to the charged Yukawa couplings may induce quantum oscillations 
amongst the lepton asymmetries. This interesting result, in combination with the temperature de- 
pendent relation between chemical potentials, makes 0(1) changes to the final value of the baryon 
asymmetry. At temperatures such that the quantum oscillations can be neglected, we find that the 
usual washout terms in the Boltzmann Equations are modified, and that there are additional terms 
proportional to flavour asymmetries. 

The upper bound on the CP asymmetry produced for each flavour is weaker than the bound 
applicable to the total CP asymmetry. In particular, it increases with the neutrino mass scale. By 
judicious choice of flavour-dependent washout rates, one can obtain a final baryon asymmetry which 
is independent of the light neutrino mass scale. Finally, given the experimental value of the Baryon 
Asymmetry of the Universe, we find a reduced lower bound on the lightest right-handed neutrino 
mass and on the reheating temperature. 
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A From a toy model to equations of motion for f l £ 

The aim of this appendix is to obtain equations of motion for a system of quantum harmonic os- 
cillators, which carry two "flavours". The oscillators are coupled to match the interactions of the 
Lagrangian in eq. ([T|). By taking expectation values of the number operator in a "thermal" state, 
this model can be reduced to a two-state quantum system. We then extrapolate the flavour structure 
of these equations of motion, to obtain "flavoured" equations of motion for the lepton asymmetry 
number operator in the early Universe. We will discuss the field theory derivation of such equations 
in a later publication [15] . 

Consider a flavour-dependent number operator, which we call /, for a system consisting of two 
simple harmonic oscillators labelled e and \x: 

(40) 



11 For complex (7^3 = U T 3 = e' s /y/2, with <j> very close to — 7r/4 and siring = — cosh r\ sin 6/(1 + cos<5)we obtain small 
Ai M . 
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where the commutation^ 2 ! relations are [a ej Oe] = 1, [a^, aji] = 1, and the hamiltonian is Hq = 
(oj e ata e + LO a a^a^)I. The number operator evolves according to the Hamiltonian equations of motion: 
^ = +i[H, /]. Since our Hamiltonian conserves particle number, we can take expectation values — for 
instance, in a thermal bath — and reduce our toy model to a two state system, described by / = (/), 
satisfying 



( 1 
dt 



(0J e - 0J u )f 



e/j 



(41) 



(w„ - LV e )f^ e 

Notice that quantum mechanically, / is the density matrix of the two-state system. If at t = the 
system is created in some state, then the probability to be found at a time t later in this state is 
Tr{/(0)/(i)}, which can oscillate. 

The excitations of the oscillators can be imagined to be the lepton asymmetry, carried by particles 
of energies uj a ~ u; e ~ E. Then uj a — tu e ~ m 2 u (T)/E, where m 2 (T) = /i 2 T 2 /16 is the contribution of 
hn to the "thermal mass" of t a in the early Universe. Suppose the production and washout of the 
lepton asymmetry are treated as initial condition and subsequent measurement on the system. That 
is, we imagine to produce the asymmetry in some linear combination of e and \i, allow it to evolve, 
and then at a later time turn on the inverse decays. Then the inverse decay rate could be reduced 
(or increased) because the asymmetry changed flavour during time evolution (analogous to a survival 
probability in neutrino oscillations). 

We have neglected the thermal masses oc AjA^ because we anticipate that its effects are negligeable 

at temperatures ~ 10 9 GeV: if \h u \ 2 3> |A| 2 , the A-contribution to Hq is small and can be ignored. If 
\h u \ 2 <C |A| 2 , then the asymmetry is produced in the flavour direction A, so a thermal mass for this 
direction will not cause oscillations (asymmetry in a mass eigenstate). For a detailed discussion of 
"neutrino flavour effects" (due to A) in resonant leptogenesis, see [13] . 

The decays and inverse decays due to the Yukawa interactions can be included by perturbing in 
an interaction Hamiltonian, Hj. We first consider the production and washout of leptons, due to A, 
in a model of four harmonic oscillators: N,H, and two leptons flavours {li}- Later we will treat the 
N and H number densities as background fields, and take the difference with respect to the equation 
for anti-lepton flavours. We are looking for behaviours analogous to the lepton number production 
and washout in the decay /inverse decay N — > Hi. 

We suppose the H and i are massless and N is massive, and ignore the free Hamiltonian (con- 
sidered previously) which does not change the particle numbers. The interaction Hamiltonian is 

Hi = AjtajyO/ra* + \%a N a} H a & \ (42) 

We perturbatively expand the Heiseberg equations of motion to get that 



Q 

—ff = -[Hi, [Hi, ff]) = -\* k \i[a N (J H af ', a) N a H a{] + \ k X f j [a} N a H a'l, a N a) H d^ 

\,Kifu.lt J - f N f H s kj - f N ft - f N P k ) 



+{f N 5 lk + f N f H S lk + f N ft ~ fHf! k )h\* . (43) 

This looks like inverse decays and decays of N, and the middle two terms of each parentheses we 
drop because they would not be allowed kinematically for N decay in the early Universe. A similar 
calculation El can be performed for the muon Yukawa coupling h, with 

Hi,h = hko^a^a'l + h* k a H a^ca k] ', (44) 



12 In studying flavour, we consider "bosonic" leptons. The generalization to the realistic case of fermionic leptons is 
straightforward. 

13 We neglect possible interference terms between Hi and Hi,h, because the expectation value of terms such as a^aN 
will vanish. 
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and using the kinematic constraint that the processes H «-> t^n are allowed, because the thermal 
mass of Higgs has contributions from the top. Taking expectation values, this gives in the charged 
lepton mass eigenstate basis: 



a_ 

Of 



fee f^/J- 

H J e 



— ffi c 

+2f H 












fee 

h 

fixe 
J£ 



h 
it. 



fee 

h 

rfie 
if. 



ne/J. 

h 

J P. 










\h I 2 



\h I 2 
I n fi | 



It is reassuring that the i-independent equilibrium solution of this equation imposes fn°fg ■ 
Combined with eq. (|4ip. this gives damped oscillations for the off-diagonals: 



(45) 

fH = 0. 



dfT/dt 



-l\Up. 



i\ h v\ 2 fv c )f, 



en 



so we can extrapolate that exchanging doublet with singlet leptons destroys the coherence between 
different doublet flavours. 

Equations for the anti-lepton density can similarly be derived. For simplicity, we take H and 
N to be their own anti-particles, so there are now eight coupled oscillators: N,H,i l ,i' l fi c and £t c . 
The expectation value of the difference between the two equations, again in the charged lepton mass 
eigenstate basis, is 



^fX = -\iXtf H f%-f H fi?\ k \* + 2(f N -f N q )eV 

-id* - Uj )fi - ^£(/ M c + fcx/^ + srfg) 

-^fu^ - u°Wi + ffy^i + Wt + ffr j ), (46) 

where k, i, j can be e or //, and i, j are not summed. The purpose of this equation is to motivate the 
flavour index structure we use in the text. We now briefly discuss the various terms. 

• The first two terms describe washout by inverse decays. By construction, there is no washout 
by AL = 2 scattering (we will include this by hand). 

• The third term, which is oc e, has been obtained artificially. It is well-known in field theory 
that no asymmetry is created at second order in A; the difference T(N — > Hi) — T(N — > Hi) 
is proportional to the imaginary parts of the loop amplitude times the loop coupling constants. 
We write the latter as (<5A)j, replace the "tree-level" Aj by a "tree + loop" vertex A^ = Aj+(<5A)j, 
so T — f becomes the matrix in flavour space [(5X)i\*— Xi(5X)*]. This gives 

en = %{Xi(8X)* - (6X)iX*j} = %{X iK * k X k - Kik X* k X*} (47) 
where = [m„]j«/ti 2 . Using the equilibrium condition 

2f N q [e} = f -f{[e][h + fj] + [fi + fj M), 

(where square brackets are matrices in flavour space), we obtain a production term 2(/jv + 
f e ^)e 13 ■ This is incomplete, because it would produce an asymmetry in thermal equilibrium. 
The resonant part of AL = 2 scattering should be subtracted [36], which ensures that Tr /a 
vanishes in thermal equilibrium. We work in the charged lepton mass eigenstate basis, and 
follow [36] to obtain eq. (06]). 
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• The remaining terms of eq. (|46p describe the effects of the muon Yukawa coupling. The 
asymmetry oscillates in flavour space, as described at the beginning of the section, because the 
off-diagonal flavour matrix elements have time-dependent phases. 

• The muon Yukawa coupling also causes the off-diagonal elements to decay away, via the second 
to last term. To get the correct relative normalisation between oscillation and decay in the 
early Universe, we will extract the rates from the real and imagainary parts of the thermal 
propagator. 

• The last term has minimal effect on the off-diagonal terms, as one can see by solving eq. (|45p . 
However, combined with the second to last term, it drives the diagonal asymmetry Yj\i^ to 
the correct relation with Yb/3-£ m > when the muon Yukawa coupling is in equilibrium o For 
simplicity, we drop this term, assuming its effects could be approximated as time-dependent 
components in the A-matrix [11] relating Y^. = A Jl Y^/ 3 _x Ji ■ 

A ij (t) = A^ >e - r ^ + ,4^(1 - e- r "*), (48) 

where Aj, y are the matrix elements before comes into equilibrium, Aj, < are the matrix 
elements afterwards, and is the interaction rate associated to the muon Yukawa coupling, 
see also appendix C. 



Added Discussion 

The effect of the charged lepton Yukawa in equation (|46p can be clarified by studying the evolution 
of i™ = — ^^-(/^c — fa), where hi is the muon Yukawa coupling in some arbitrary basis for 
lepton doublet space. The trace of F is the total lepton number stored in ex,, hl and /j,r, so studying 
F in the toy model is analogous to replacing the asymmetry in lepton doublets with the asymmetry 
in B/3 — L a [11] in usual leptogenesis calculations. Neglecting the first four terms of eqn (06]), which 
describe the interactions of ./V and the oscillations due to the muon Yukawa, the last two terms of 
eqn (|46p can be written 

| [/a] = ~{[hh% [/ A ]} (/ * C + f ^ - {[ffi], [f £ + ff]} (/ *° " M . (49) 

where h is a colomn vector, and objects in square brackets are matrices in doublet space. The rate 
of change of the /x c asymmetry = f^c — f-^c is 

= -\h\ 2 (U° + fr)f!?-\h\ 2 fAr(Ji + fi) Mi 

So the equation for the lepton asymmetry i™ , in the flavour basis, is 

= - ({m, if A ]} - 2 w • [/ a ] • k) (/ ^ c + /?zc) 
- ({m, \h + fj ]} - 2 w . [ife + f T )} . ^ (^i_m . (50) 

where the second line is of the same order as the first (see eqn (|46p ).As discussed in the last point 
after eqn (|46p . an approximate equation for F 33 can be obtained by expressing /a^ c and in terms 



14 In our toy model, where H is a real scalar, this means (/ M c — f^c) = f^ 1 . In the SM, there is Higgs asymmetry 

rmrni. 
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of F n using a time-dependent A-matrix. To obtain equation (jlip . we simply drop the second line. 
In the flavour basis, eqn (|50p can be written as 



0_ 

Of 



fee f e / 1 
J A J A 

f lj,e f 

J A J A ~ J^ c 



o ft 

fA 



\h I 2 



/r + /r ° 



As expected, the Yukawa coupling h causes the flavour off-diagonal asymmetries to decay away, as 
it did for the off-diagonal number densities in eqn (j45f) . It has no effect on the asymmetry in either 
flavour. If the matrices F and \hiv] = H a a a are expanded on Pauli matrices, as in equation (|18p . 
this behaviour can be modelled as a triple cross product: H x H x F. A single cross product, 
as in the orginal version of the paper, obviously does not cause the off-diagonals to be damped as 
dF ± /8t oc \h\ 2 F ± . 



B AL = 2 terms 



Following the same procedure of section [4] we can include the contributions arising from the AL = 2 
exchange mediated by RH neutrinos. In matrix form, the additional terms on the RHS of eq. ()23[) 
arising from the s-channel contribution would be 



RHS of eq. {23 




eeee 
7AL=2 



Ufiee 
IAL=2 



Tree 
7AL=2 



yec 



V TT 



eenn 
7AL=2 



eerr 
7AL=2 



'AL=2 



/AL=2 



TT[1U 

'AL=2 



TTTT 

7al=2 



Vee 
L 



Vee 



VTT 

L 



V TT 



Tt[7aT= 2 ^ J1 



(51) 



where 

itl 2 = e -w. = e i^n^f^- (52) 

The sum is over the RH exchanged neutrino eigenstates and the function Fjj is the kinematical factor 
corrected according to [5] . We define the matrices 



[7. 



AL=2J 



..nee 
7AL=2 



UUfl 

>AL=2 



ryllTT 

!AL=2 



(53) 



and 



[7al= 2 ] 



>AL=2 



<AL=2 



>AL=2 



(54) 
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Each of which can be decomposed in terms of the matrices /, A3, A§: 



[llL=2\ = \lAL=2 I + 7^AL=2,3A 3 + 7AL=2,8 A 8 



(55) 



and similarly for 7al = 9 • The interference terms can also be dealt with in terms of the quantities 



(56) 



KJ 



where Gjj is a kinematical factor of the same kind as Fjj [5] . Thus taking the trace the contribution 
to eq. ([2j 



RHS 



[7AL=2 + 7AL=2 + 7AL=2 + 7AL=2 + 7al=2 + 7al= 2 ] y l 



1 

+ 3 
+2 

+2 

+2 

+2 



■ e(IJ) u(IJ) r(IJ) e(IJ) MU) r(IJ) 

7AL=2 + 7AL=2 + 7AL=2 + 7AL=2 + 7al=2 + 7al=2 



Yl 



7AL=2,3 + 7AL=2,3 + 7AL=2,3 + 7aL=2,3 + 7aL=2,3 + 7aL=2,3 

e(IJ) u(IJ) + t(IJ) e(IJ) MI J) , ,r(JJ) ' 

7AL=2,3 + 7AL=2,3 + 7AL=2,3 + 7aL=2,3 + 7aL=2,3 + 7aL=2,3 



7AL=2,8 + 7AL=2,8 + 7AL=2,8 + 7aL=2,8 + 7aX=2,8 + 7aL= 



:2.8 



e(IJ) u(JJ) r(JJ) e(JJ) ,m(/J) , ,t(IJ) 

7AL=2,8 + 7AL=2,8 + 7AL=2,8 + 7aL=2,8 + 7aL=2,8 + 7aL=2,; 



*3 

^8 
Is 



(57) 



We next show that for low temperatures (z > 1), the AL = 2 terms can be neglected with respect 
to the inverse decays. To simplify the discussion we focus on the ID and AL = 2 contributions to a 
single flavour, say elements 11 of the correspondings matrices that enter the BE. For z>l 



lD Y L r E/ 9 ■,/9 M l X le*le Y L 



IT 



5/2 z 3/2- 



(58) 



and the AL = 2 terms give, 



RHS of Eq.m oc ( 2|A le | 2 | ( £ A lfe | 2 ) Y£ e + |A le | 2 |A lM | 2 |y L w + |A le | 2 |A lr | 2 |Y L " 



(59) 



This means that for temperatures z 55 z c , we can neglect the contribution to the BE from the 
lepton number violating exchange of the RH neutrinos as the comparative strength of the inverse 
decays suppression is ~ e~ z while the suppression from the AL = 2 is 



Tr 5 / 2 Z 9/2 



1 J- T 



\k\ + \^lfj.\ 2 T?77 + I^1t| 2 



vee 
1 L 



(60) 



This implies that for |Aife| 2 ~ 1 only for z ^ 15 will the AL = 2 terms dominate over 773. For 
smaller values of Ai/%, as constrained by the departure from thermal equilibrium and neutrino masses, 
the AL = 2 terms are always neglible for the more interesting values of z. That is, if the freeze out 
temperature is denoted by Zf, then for the regions in parameter space for which Zf ^ z c the produced 
asymmetry is not affected by the AL = 2 washout terms and thus can be neglected. 
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C Chemical Equilibrium 



An important point which was pointed out in ref. is the relationship between the lepton asym- 
metry for each flavour Li and the B/3 — Li = Aj asymmetry, which is what is effectively conserved by 
the sphaleron interactions, from chemical equilibrium equations. More recently ref. [37] has also dis- 
cussed issues related to this. For completeness we include the corresponding matrices which relate Aj 
to the chemical potential of the left-handed leptons, [A Ui , depending on which of the charged Yukawa 
couplings are in equilibrium as discussed in the main text of the paper. Note that the relationship 
between \i Vi and Li also varies depending on the corresponding charged Yukawa couplings that are 
in equilibrium. 

When T T 3> H, and sphaleron interactions are also in equilibrium, then 




_22 _4 _4 
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9 9 9 

4,5 1 4,5 1 4,8 1 

9 "r" 3 3JV+3 9 ~r 3 3N+3 9 3 3N+2 







( ^ \ 




1 
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(61) 



where N is the number of generations. For the range of temperatures we have focused on in the main 
part the relation between these quantities is modified to: 




22 4 4 
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4 I r 1 4 o i 8 1 4,8 1 II J (f , 2) 

9^3 3iV+4 9 ° ^ 3 3JV+4 9^3 3JV+4 I I AS. I ■ l uz J 





( ^ \ 
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9^3 3N+4 9^3 3AT+4 9 u ^ 3 3A^+4 

Finally, when all charged Yukawa couplings are in equilibrium, the relation is simply 

Ai =B + 54A^)S^-^,. (63) 

j 
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